Source of Acquisition 
NASA Goddard Space Flight Center 


What does reflection from cloud sides tell us about vertical distribution of cloud 

droplet sizes? 

Alexander Marshak 1 , J. Vanderlei Martins 1,2 , Victor Zubko 3 , Yoram J. Kaufman 1 

'NASA/GSFC, Climate and Radiation Branch, Greenbelt, MD 
2 Joint Center for Earth Systems Technology, UMBC, MD 
3 Advanced Computation Technologies, Inc., Lanham, MD 

Public Summary 

The abundance and nature of the aerosol particles determine the characteristics of the 
distribution of cloud droplets condensing in the cloud. These, in turn, influence cloud structure, 
freezing level and the height the cloud needs to reach for rain to be started. Thus, in order to 
understand the effect of aerosol on cloud development and precipitation, we need to measure the 
vertical profile of the droplet sizes and determine the level of freezing in the cloud. Since cloud 
development varies for different meteorological and geographical conditions, a satellite mission is 
the only way for getting global information. 

The recently proposed CLAIM-3D satellite mission (cloud aerosol interaction mission in 
3D) suggests to measure vertical distribution of cloud droplets by retrieving them from the solar 
and infrared radiation reflected or emitted from cloud sides. This paper discusses how to interpret 
satellite measurements from the cloud sides assuming that we measure the reflected sunlight from 
the cloud sides and top at two wavelengths: one where solar radiation is not absorbed and one with 
liquid water absorption of solar radiation. 

Scanning the cloud sides is fundamentally a three-dimensional (3D) problem. Blindly 
applying traditional cloud retrieval methods that do not take into account the cloud horizontal 
inhomogeneity and assume that the satellite always sees the cloud top, may produce erroneous 
results that could lead to misinterpretation of the physics involved in cloud development. Here, 3D 
radiative transfer is used for interpreting the observed cloud reflectances. As a proof of concept, 
the paper shows a few examples of radiation reflected from cloud fields generated by a simple 
stochastic cloud model with the prescribed vertically resolved microphysics. To retrieve the 
information about droplet sizes, the probability density function of the droplet size distribution 
rather than a fixed value of effective droplet radius is used. 
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Abstract 


Cloud development, the onset of precipitation and the effect of aerosol on clouds depend on 
the structure of the cloud profiles of droplet size and phase. Aircraft measurements of cloud 
profiles are limited in their temporal and spatial extent. Satellites were used to observe cloud tops 
not cloud profiles with vertical profiles of precipitation-sized droplets anticipated from CloudSat. 
The recently proposed CLAIM-3D satellite mission (cloud aerosol interaction mission in 3D) 
suggests to measure profiles of cloud microphysical properties by retrieving them from the solar 
and infrared radiation reflected or emitted from cloud sides. 

Inversion of measurements from the cloud sides requires rigorous understanding of the 3- 
dimentional (3D) properties of clouds. Here we discuss the reflected sunlight from the cloud sides 
and top at two wavelengths: one nonabsorbing to solar radiation (0.67 (im) and one with liquid 
water efficient absorption of solar radiation (2.1 jam). In contrast to the plane-parallel 
approximation, a conventional approach to all current operational retrievals, 3D radiative transfer 
is used for interpreting the observed reflectances. General properties of the radiation reflected 
from the sides of an isolated cloud are discussed. As a proof of concept, the paper shows a few 
examples of radiation reflected from cloud fields generated by a simple stochastic cloud model 
with the prescribed vertically resolved microphysics. To retrieve the information about droplet 
sizes, we propose to use the probability density function of the droplet size distribution and its first 
two moments instead of the assumption about fixed values of the droplet effective radius. The 
retrieval algorithm is based on the Bayesian theorem that combines prior information about cloud 
structure and microphysics with radiative transfer calculations. 
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Abstract 
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dimentional (3D) properties of clouds. Here we discuss the reflected sunlight from the cloud sides 
and top at two wavelengths: one nonabsorbing to solar radiation (0.67 pm) and one with liquid 
water efficient absorption of solar radiation (2.1 pm). In contrast to the plane-parallel 
approximation, a conventional approach to all current operational retrievals, 3D radiative transfer 
is used for interpreting the observed reflectances. General properties of the radiation reflected 
from the sides of an isolated cloud are discussed. As a proof of concept, the paper shows a few 
examples of radiation reflected from cloud fields generated by a simple stochastic cloud model 
with the prescribed vertically resolved microphysics. To retrieve the information about droplet 
sizes, we propose to use the probability density function of the droplet size distribution and its first 
two moments instead of the assumption about fixed values of the droplet effective radius. The 
retrieval algorithm is based on the Bayesian theorem that combines prior information about cloud 
structure and microphysics with radiative transfer calculations. 
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1. Introduction 


Investigation of cloud development and the onset of precipitation are essential to 
understand the role of clouds in the hydrological cycle and the effect of pollutants on clouds and 
precipitation (Ramanathan et al., 2001). It also advances our understanding of the feedback of 
clouds on climate change and the aerosol indirect forcing of climate through cloud modification. 
Therefore, we have to resolve the vertical distribution of cloud droplet sizes and determine the 
temperature of glaciation for clean and polluted clouds (Andreae et al., 2004). Knowledge of the 
droplet vertical profile is also essential for understanding precipitation (Rosenfeld, 2000, 
Rosenfeld and Ulbrich, 2000). In an accompanied paper, Martins et al. (2006) suggest a satellite 
mission to derive profiles of the cloud microphysics using observations of the cloud sides. Here 
we show a methodology, based on 3-dimensional (3D) cloud properties to retrieve the cloud 
profiles from the new satellite measurements. 

So far, all existing satellites either measure cloud microphysics only at cloud top (e.g., 
Moderate Resolution Imaging Spectrometer (MODIS), see Platnick et al., 2003) or give a vertical 
profile of precipitation sized droplets (e.g., CloudSat, see Stephens et al., 2002). Note that the 
combination of millimeter-wave radar reflectivity measured by CloudSat with MODIS (on Aqua) 
measurements of solar radiance will be able to provide cloud droplet size vertical profiles but 
under some strong assumptions of given number concentration and droplet size distribution. 

Except for Polarization and Directionality of the Earth’s Reflectance (POLDER) that 
retrieves cloud droplet effective radius at the very top cloud layer (with an optical thickness of 1) 
from polarization measurements of the reflected light (e.g., Breon and Golub, 1998, Breon and 
Doutriaux-Boucher, 2005), all operational retrievals of cloud droplet size are based on spectral 
observations (e.g., Nakajima and King, 1990). For MODIS, cloud optical thickness, t, and droplet 
effective radius, r e , are simultaneously derived from various two band combinations: typically one 
water-absorbing band {1.6, 2.1, or 3.7 jum} and one nonabsorbing band {0.65, 0.86, or 1.2 pm} 
(Platnick et al., 2003). Since water absorbs differently in the three MODIS absorbing bands, the 
less absorbing 1.6-pm band and the more absorbing 3.7-pm band complement to the 2.1 -pm band 
for assessing the vertical variation of r e in the upper portion of the cloud (Platnick, 2000; Chang 
and Li, 2002). However, these variations are not sufficient to resolve the vertical distribution of 
cloud droplet sizes from cloud base to cloud top. What is if one would measure the vertical 
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profiles of the cloud microphysical properties by retrieving them from the solar (and infrared) 
radiation reflected (or emitted) directly from cloud sides? 

Note that all existing operational retrieval algorithms are based on the plane-parallel 
approximation that does not take into account the cloud horizontal inhomogeneity. In terms of 
cloud aspect ratio, A=L/h (where L and h are horizontal and vertical dimensions of a cloud, 
respectively), the main plane-parallel assumption used for any remote sensing retrieval is that A is 
infinitely large and that the satellite always sees the cloud top. Hence, a pair of reflectances at the 
nonabsorbing and absorbing bands indicates how optically thick (thus estimates x) and how 
absorbing (thus estimates r e ) clouds are (Nakajima and King, 1 990). 

It is well understood that finite isolated clouds of various shapes and sizes can have 
absolutely different radiative properties than their plane-parallel counterparts. Davies (1978) 
represented an isolated cloud as a cuboid of given dimensions. In this case, the incident solar 
beam hits not only the top of the cloud but also one or two of its sides. As an alternative to the 
plane-parallel model to simulate cumulus clouds, recently Davis (2002) used a spherical turbid 
medium. For his spherical cloud, he was able to derive analytically the transmitted and reflected 
fluxes in terms of the cloud optical diameter. He showed that these results could be used to 
estimate the cloud optical diameter from radiances reflected from dark and bright sides of clouds. 

In general, if one releases the assumption that the aspect ratio A is infinitely large then, in 
addition to cloud tops, a satellite-based observer will likely see cloud sides. Because of a variety 
of possible aspect ratios and cloud geometrical shapes, the situation seems to be out of control and 
measured data camiot be correctly interpreted in the sense of cloud properties. Similar to the 
plane-parallel approximation, in order to bring the retrieval back under control we have to make 
simplifying assumptions. The main assumption for cloud side remote sensing is that regardless of 
the aspect ratio, cloud geometrical shape and its microphysical structure, a pair of reflectances at 
nonabsorbing and absorbing bands determines a distribution of droplet sizes. Note that this is an 
assumption rather than a statement since it can’t be checked with the model calculation and 
inversion for all cloud types. Also note that here we are talking about the distribution of droplet 
sizes (with mean and standard deviation) rather than a single value. Finally, together with the 
brightness temperature this assumption allows us to estimate a vertical profile of droplet (particle) 
sizes (Martins et ah, 2006). 
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Of course, the above assumption will not work for all cloud types like the plane-parallel 
approximation does not work for all clouds either. Here we will consider only optically thick 
clouds (t> 40) with relatively small aspect ratio (L/h < 2-5). We will further make some 
additional limitations regarding the satellite viewing angles. In order to see a sufficient amount of 
cloud sides, the viewing zenith angles, 0, will be limited to more oblique angles of 0 > 45°. For 
simplicity here we will be considering only “backward” directions, i.e., cp~(p 0 where cp 0 and 9 are 
solar and viewing azimuthal angles, respectively. Under these rather strong limitations, the paper 
proves the concept of a possible retrieval of the distributions of droplet vertical profiles using three 
bands: nonabsorbing, water absorbing and brightness temperature. The latter is associated with the 
measured height and is discussed in the companion paper (Martins et al., 2006). 

The plan of the paper is as follows. Section 2 briefly discusses the main radiative transfer 
features of the reflectance from cloud sides based on a single homogeneous cloud. To generalize 
these results to a more realistic horizontally inhomogeneous cloud field, Section 3 describes simple 
stochastic and microphysical models used to simulate a variety of cloud fields. With the help of 
two wavelengths at 0.67 and 2. 1 pm, Section 4 demonstrates the retrievals of the distribution of 
droplet sizes from the measurements of radiation reflected from the cloud fields simulated in 
Section. 3 . At the end of Section 4, this approach is generalized in the terms of Bayesian retrievals 
(McFarlane et ah, 2002, Evans et al., 2002). Finally section 5 provides general discussion and 
summarizes the results. 

2. Radiative transfer calculations 
2.1 3D radiative transfer tools 

There are two 3D Radiative Transfer (RT) tools that dominate atmospheric radiation 
applications and are currently the only available options for solving complex RT problems: the 
Spherical Harmonic Discrete Ordinate Method (SHDOM) of Evans (1998) and the Monte Carlo 
(MC) method (Marchuk et al., 1980). When many radiative quantities are required, e.g., the 
radiance field across cloud top, SHDOM is much faster than MC, but its errors (and limitations) 
are harder to interpreter, especially for optically thick and highly variable media around cloud 
edges. Since the rule-of-thumb in using SHDOM requires the optical path across a grid cell to be 
of order of one, its solution may be not accurate for horizontally and vertically thick clouds. 
Moreover, SHDOM (tri)linearly interpolates the extinction between grid points; thus it may have 
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some problems when reflectance from cloud sides of optically thick clouds is calculated. Anyway, 
in this study we used both MC and SHDOM; for several key calculations both methods were 
applied simultaneously to the same set of cloud parameters to intercompare and validate the 
results. To the best of our knowledge, the results shown in the paper are numerically accurate. 
(For the detailed description of both the SHDOM and MC methods, see Evans and Marshak, 
2005). 

2.2 Main radiative transfer features of the reflectance from cloud sides 

Using a 3D Monte Carlo radiative transfer code (Evans and Marshak, 2005), we calculated 
reflectance at 0.67 and 2.1 pm wavelengths from a single homogeneous cloud. The cloud top was 
simulated by an infinitely long rectangular with width L= 12.8 km and height h= 2 km (Fig. la). 
Cloud vertical optical thickness, x, varied from 20 to 160 and droplet effective radius was assumed 
a constant r e =T0 pm. The cloud was illuminated at solar zenith angle (SZA) 9o=60° along cloud 
inhomogeneity and observed at viewing zenith angles (VZA) 0 = 45°-70° from the illuminated side 
of the cloud. The reflectances are plotted in panels (b) and (c) on Fig. 1. The horizontal axis 
shows the distance to the cloud edge outside the cloud (negative x-values) and the distance from 
the cloud edge inside the cloud (positive x-values). The cloud edge is located at x=0. For 
example, a cloud side, viewed at 0=70°, can be seen /zxtan(0)=5.5 km away from the cloud 
(negative 5.5 km). Thus negative x-values correspond to radiation reflected from a cloud side 
while positive x-values to radiation reflected from a cloud top. Here are the main features that can 
be observed from these two panels. 

• Reflectance from a cloud side at 2.1 pm is saturated starting from t= 40 while reflectance at 
0.67 pm does not reach the level of saturation at all or will be saturated only at very large 
values of cloud optical thickness t. The maximum 2.1 pm reflectance from cloud sides, 
7 S ide(9o,9), depends on both SZA, 0o, and VZA, 0. It can be estimated as 

/side(Oo, 0) < /p P (9O°-0o, 9O°-0, (p-q>o=l 80°, x=-) cos(9O°-0 o ) / cos(0 o ), (la) 

where / pp (9o,0,(p"(port) is the cloud top reflectance calculated using the plane-parallel 
approximation (Stamnes et ah, 1988). For example, for 0o=9=6O°, the 2.1 pm reflectance 
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4ide(60°, 60°) < /p p (30°, 30°, cp— (p 0 — 180 °, t:= 1 60) cos(30°) / cos(60°) = 0.782, (lb) 
as seen in panel b. 

• The more oblique viewing zenith angle© (or the larger cloud side, h) the wider the area of 
maximal reflectance at 2.1 pm (panel c). 

• For optically and geometrically thick clouds, the reflectance from cloud side near cloud top is 
smaller than the one reflected from the middle of the cloud side. This effect is much more 
pronounced for 0.67 jam than for 2.1 pm. 

• For thick clouds, starting from a few optical depths away from cloud edges, reflectance from 
cloud top at 2. 1 pm is well approximated by the plane-parallel approximation. Depending on 
the extinction coefficient, it is not always the case for reflectance at 0.67 pm. At both 
wavelengths reflectance from cloud top increases towards the illuminated side and decreases 
towards the shadowed side. 

• Finally, the number of measurements from cloud side is equal to /i><tan(0)A where s is the 
horizontal resolution of a radiometer. For example, if h=2 km, 0=70°, and 5=0.1 km, there will 
be 55 cloud side measurements. 

All of the above radiative transfer features will be observed by analyzing the reflectance from 
more complex cloud fields. 

2.3 Reflectance from cloud sides for a cloud with variable droplet sizes 

Figure 2 shows an example of reflectances from cloud side and cloud top for the same two 
wavelengths (0.67 pm and 2.1 pm) but with droplet effective radius r e increasing linearly with 
height from 5 pm (at the cloud base) to 25 pm (at the cloud top). Cloud geometrical thickness h = 4 
km and cloud optical thickness is i=80 (thus extinction coefficient is 20 km' 1 ). With horizontal 
resolution 5=0.25 km and VZA 0=45, there are Man(0)/5=16 cloud side “measurements.” As for a 
simple example in Fig. 1, /o .67 reaches its maximum near cloud top (actually about 1 km from the 
cloud top) where yet most of the photons are reflected back from the cloud side without either 
transmitting through cloud and escaping from cloud base or reflecting from cloud top. Unlike in 
the previous example, the horizontal size L of a cloud is only 6.5 km and with the extinction 
coefficient 20 km" 1 this is not sufficient to reach a stable plane-parallel regime at cloud top. As a 
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result, io .67 keeps decreasing from the illuminated cloud edge to the shadowed one. In contrast, 
I 2 . 10 has a flat plateau of 5 km across where the 3D reflectance perfectly matches the plane-parallel 
one. Because of increasing droplet sizes with height, the maximum is reached much lower than in 
case of conservative scattering. It is around 1 km from cloud base where r e =9-11 pm. With 
farther increase of r e , reflectance 72 . 10 drops fast and reaches a flat plane-parallel level already at 
the cloud top (r e = 25 (im) about 1 km from the cloud edge. 

The study of reflectance from an isolated finite-size cloud is not new and has begun yet in 
early 70-ies (see, e.g., McKee and Cox, 1974, Davies 1978 and 1984). As it is seen from Figs. 1 
and 2, cloud side reflectances at the two (water-absorbing and nonabsorbing) wavelengths, have 
well-determined features. Not unlike their cloud top counterparts in the plane-parallel 
approximation (Nakajima and King, 1990), the combination of these two reflectances can be 
mapped into retrievals of cloud optical (t) and microphysical (r e ) structure. The key question here 
is whether these features survive if applied to realistic cloud fields rather than a single isolated 
horizontally homogeneous cloud. Next we briefly discuss cloud models used in this study. 

3. Cloud models 

Realistic 3D cloud fields, as an input in radiative transfer calculations, can be obtained 
from either dynamical or stochastic cloud models. For the purpose of this paper (to learn what 
reflection from cloud sides tells us about vertical distribution of cloud particles), a choice of model 
is not very crucial. The main requirements for a model were set as to have a field of several joined 
and disjoined clouds with the prescribed (observed) mean, standard deviation and correlation 
function of variable cloud optical thickness x(x,y) with a desired cloud fraction A c and cloud top 
height h(pc,y). Having some experience in stochastic cloud modeling (e.g., Marshak et al., 1994, 
Prigarin and Marshak, 2005), we selected a broken cloud version (Marshak et al., 1998) of a 
simple fractionally integrated cascade model (Schertzer and Lovejoy, 1987) that generates cloud 
fields with a given power spectral exponent, mean and standard deviation of cloud optical 
thickness. To correlate x(x,y) with h(x,y), we generated independently a t(x,y)-field and the mean 
photon free path field l(x,y). The cloud geometrical thickness field (assuming cloud base to be a 
constant) is thus a product between the optical depth and the mean free path fields, 

h(x,y) = x (x,y) * l(x,y) (2) 
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Figure 3 illustrates one realization of a cloud with an array of optical and geometrical thicknesses. 
Though it might not look very realistic, it preserves the observed correlation function in both 
optical and geometrical thicknesses. 

After cloud structure, cloud microphysics is perhaps the most important cloud model 
feature needed for radiative transfer calculations. For simplicity and for more straightforward 
interpretation of the simulated radiative transfer results, we made two assumptions: 

• cloud droplets grow linearly with z, i.e., 

r e (z‘, x,y ) = a(z-Zo) + b,z 0 <z< h(x,y), (3a) 

• the extinction coefficient <7 ex t does not depend on z, i.e., 

o'ext (z-,x,y) = a ext (x,y). (3b) 

Note that under some general assumptions (e.g., Platnick, 2000), cloud liquid water content (LWC) 
is proportional to a product of the density of liquid water, p, cube of the droplet effective radius, r e , 
and the total number of droplets in unit volume, A, 

4 

LWC(z;x,y)~-7tpr 3 '(z;x,y)N(z\x,y). (4a) 


Cloud LWC is also related to x, r e , and p as (Stephens, 1994, p. 219) 




x{x ' y)= Y P I 


LWC(z;x,y ) 

cfe Ay) 


dz. 


(4b) 


Therefore, with the assumptions (3a)-(3b), A changes with vertical coordinate z as r e " 2 , namely, 

<x,y) 1 1 cr ext (x,y) 


N(z;x,y)-- 


At the cloud base for z=zq, we get 


N(z 0 ;x,y )■ 


2nh(x,y) r 2 (z;x,y) 2 n r 2 (z;x,y ) 


l ^(Ay) 


(5a) 


(5b) 


2k b 2 

For example, for r c (z 0 ) = b = 5 pm and a ext = 20 km 1 one gets N(z, 0 ) = 127 cm" 3 . If at the cloud top 
r e = 25 pm then Eq. (5a) yields N(z 0 ) = 5 cm" 3 . Figure 4 shows an example of vertical profiles for 
cloud liquid water, LWC (in g/m 3 ), total number of drops, A (in cm" 3 ), effective radius, r e (in pm), 
and extinction coefficient, G ext (in km" 1 ). We see that, for each x and y, with z increasing from 
cloud base zo to cloud top h(x,y), LWC and r e increase linearly, A decreases as z 2 , and o ext is 
constant. 
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4. Proof of concept 

Figure 5 shows an example of reflectances from a 16 by 16 km cloud field illuminated at 
0o=6O° [from South (bottom of the image)] and viewed at different viewing angles: 0= 0°, 20°, 45°, 
and 70° (also from South). The cloud is 4 km thick; for illustrative purposes, the cloud top is flat. 
Droplet effective radius grows linearly with height from 5 to 25 pm; thus in Eq. (3a), a = 5 and b = 5 
pm. 

The two upper plots show nadir angle observations. As illustrated in Figs. 1 and 2, we see 
that at 0.67 pm, the cloud tops at the illuminated cloud edges are much brighter, whereas the cloud 
tops at the opposite ends look darker then in the rest of the area. At 2.1 |im, cloud tops are 
homogeneous except may be the first 0.5 km away from the illuminated cloud edges. With 
increasing viewing angles, we start seeing illuminated cloud sides that are brighter than their cloud 
top counterparts. As a result, even visually one can distinguish between cloud sides and cloud 
tops, especially at low viewing angles. Similar to Fig. 2, at 0.67 pm the reflectance from cloud 
sides reaches its maximum in the middle of the cloud while at 2.1 pm the reflectance from cloud 
sides gradually decreases starting from about 0.5-1 km (10-20 optical depths) from the cloud base. 
This decrease is a clear signature of droplet sizes that are small (5 pm) at the bottom and are large 
(25 pm) and highly absorptive at the top. 

A scatter plot on Fig. 6a is a Nakajima-King (1990) type diagram that relates cloud 
reflectances at 2.1 and 0.67 pm. The plot is based on 20 cloud fields generated as realizations of 
the stochastic cloud model described in section 3. In contrast to a traditional Nakajima-King 
scatter plot that shows only the cloud-top reflectance, most of the points on Fig. 6a correspond to 
the reflectance from cloud sides. Indeed, panels (b) and (c) illustrate the break down of panel (a) 
into reflectance from cloud sides and cloud tops, respectively. Panel (b) is much brighter than 
panel (c), i.e., much more (observed) photons have been reflected from cloud sides than from 
cloud tops. We also see from panel (c) that, since cloud droplet (particle) size increases linearly 
with height (see, Eq. (3a)), only those cloud tops that have the largest r e = 25 pm (blue dots) have 
substantially contributed to the total reflectance. Because of low VZA (0=70°), other cloud tops 
are in shadow and are barely seen by the observer. As explained in Section 2.3, at 2.1 pm the 
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cloud-top reflectances (blue dots) are the smallest. At 0.67 jim, the cloud-top reflectances have a 
wide range of values; the latter corresponds to the variety of cloud optical thicknesses as follows 
directly from the Nakajima-King (1990) theory. 

Let us now fix the 0.67 (im reflectances at four different levels (dash lines in Fig. 6) and 
build histograms of r e for different values of the 2. 1 jam reflectances. Figure 7 illustrates them 
(with a window of ±0.03 for /o.67 and ±0.02 for J 2 .i). As expected, for 7o.67-0.83, most observed 
radiances are reflected from the cloud top with r e = 25 jim. Increasing /o.67, we observe more and 
more radiances reflected from the cloud side. For the large enough 0.67 pm reflectances, / 2 .i 
saturates and, similar to the plane-parallel approximation, the retrieved values of r e become 
insensitive to the values of Io.6i- Because of the statistical nature of our retrievals, instead of a 
single value of r e , we retrieve its (conditional) probability density, p(r^\h.ei,h.\)- The mean r e can 
be calculated as 

oo 

( r e) = J r e p(.4l 0A1 ,I 2 .i)dr e (6a) 

0 

and its standard deviation a as 

<7 = Jj 0c - {r e )f p{r e \lo. 61 J 2 .i)dr e ■ (6b) 

For example, for To. 67 - 1.22±0.03 and h. i = 0.58±0.02, the mean retrieved value <r e >= 12 pm with 
standard deviation 0=2 pm. 

If, in addition to the measurements at 0.67 and 2.1 pm, one also measures the cloud side 
brightness temperature, say at 11.6 pm, each retrieved distribution of effective radius can be 
directly related to cloud side brightness temperature, thus assessing its altitude (for details see the 
companion paper, Martins et al., 2006). In other words, a combination of measurements at these 
three wavelengths can resolve the vertical distribution of cloud droplet sizes near cloud side. The 
extension of the retrieved profiles from cloud sides to the whole cloud requires an additional 
assumption of mild fluctuations of droplet effective radii along a horizontal plane at the same 
altitude z inside clouds. As discussed in Martins et al. (2006), studies of in situ measurements in 
Cumulus clouds (e.g., Blyth and Latham, 1991) and cloud models (Zev Levin, private 
communications) confirm that this assumption does not look unrealistic either. 
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Generally speaking, to retrieve a vertical profile of droplet effective radius, the above 
approach suggests using a database of stochastic cloud models and corresponding radiative transfer 
calculations of cloud reflectances at 0.67, 2.1 and 1 1.6 jim. This is similar to a Bayesian retrieval 
algorithm (e.g., McFarlane et ah, 2002; Evans et al., 2002) that combines prior information about 
cloud structure and microphysics with radiative transfer calculations, 




T ) (^0.67’^2.1’Al.6 x )p( x ) 
j P( I 0.67’ I 2.lJu.6\ X )P( X ) dX 


( 7 ) 


Here the vector x consist of cloud parameters (with ;- e ) that affect the cloud reflectances: /o. 67 , h.i 
and in 6- Function p{Io.6i,h.\,h\.e\ x ) is the conditional probability density function given vector x. 
It is directly related to our pre-calculated database — the radiative transfer simulations of cloud 
reflectances for the cloud structure defined by x. Figure 8 shows an example of the conditional 
probability function of two reflectances io .67 and / 2.1 for r e from 6 to 7 |im and from 14 to 15 |im, 
respectively. Other parameters of cloud structure (vector x) that affect calculations of io .67 and h.\ 
are described in Section 3 and in the caption to Fig. 6. Note that the distribution functions in Fig. 8 
are not necessarily Gaussian. Function p(x) is the probability density function of cloud structure x. 
In other words, based on the prior information it defines the weights to be accounted for 
characterizing the frequency of state x. The integral in the denominator on the right hand side of 
(7) is just a normalizing factor. Finally, the left hand side of (7) is the (posterior) probability 
density of having cloud structure x (including r e ) giving measurements of io. 67 , h.\ and in. It is 
related to histograms shown in Fig. 7. For details on a Bayesian retrieval algorithm applied to 
microwave radiometer and submillimeter- wave cloud ice radiometer see the excellent descriptions 
given in McFarlane et al. (2002) and Evans et al. (2002), respectively. 


5. Summary and discussion 

Knowledge on vertical distribution of droplet sizes is essential for understanding not only 
cloud development and precipitation but also the interactions between clouds and aerosols. 
Recently Andreae et al. (2004) using in situ aircraft measurements showed a strong dependence of 
various cloud properties (including droplet sizes), as a function of height in the cloud, on 
abundance of aerosol particles. How can one obtain this information globally from satellite remote 
sensing? 
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For this purpose, a new satellite mission, called CLAIM-3D (stands for ‘3D cloud aerosol 
interaction mission’) has been recently proposed (Martins et al, 2006). The CLAIM-3D mission is 
designed to advance our understanding of cloud and precipitation development by measuring 
vertically resolved cloud parameters. It is proposed to have a unique combination of extended 
wavelength range (0.38 - 12 pm), polarization, and multi-angle 3D geometry combining the best 
features from POLDER (polarization), MISR (multi-angle), and MODIS (multi-channel) to 
characterize aerosols and cloud microphysics. This paper is the first step towards possible 
interpretation of CLAIM-3D measurements of reflected from cloud sides solar radiation. 

Over the last two decades, considerable efforts have been dedicated to optical remote 
sensing of cloud properties. Using plane-parallel radiative transfer theory (e.g. Nakajima and 
King, 1990), measured radiances have been turned into science products, such as cloud optical 
depth and effective radius. If this approach is acceptable for a stratiform type clouds, it is suspect 
for clouds that are far from horizontally homogeneous (e.g., Varnai and Marshak, 2001; Iwabuchi 
and Hayasaka, 2002; Davis, 2002), especially for the clouds with a relatively small aspect ratio 
(the ratio of horizontal to vertical cloud dimensions) and well-developed cloud sides. These are 
the clouds the CLAIM-3D mission is directed for. In other words, here we target cloud side 
passive remote sensing rather than traditional cloud top remote sensing. 

However, in order to interpret the cloud side measurements, a new 3D-based cloud retrieval 
scheme should be developed. Advances in 3D radiative transfer algorithms, improved 
understanding of 3D cloud structure (Marshak and Davis, 2005), and increases in computing 
power make the time now ripe for 3D cloud retrieval. 

This paper studies the properties of radiation reflected from cloud sides at two 
wavelengths: one nonabsorbing (0.67 Jim) and one water-absorbing (2.1 pm). As a proof of 
concept, it shows that (under some general assumptions and limitations) using Bayesian approach 
(e.g., Evans, 2002) simultaneous measurements of radiances at these two wavelengths can be 
mapped into a distribution of cloud droplet sizes. Not unlike the famous Nakajima-King (1990) 
diagram that maps cloud top reflections into a pair of cloud optical depth and effective radius, a 
new algorithm based on cloud stochastic models is capable in interpreting cloud side reflections at 
0.67 and 2.1 pm in terms of cloud droplet size distribution. If the infonnation on cloud side 
brightness temperature is available, droplet size distributions can be vertically resolved. 


13 



Of course, knowledge of reflectance from the pixels surrounding each target pixel as well 
as reflectance at multiple angles will improve our retrieval making the width of the retrieved 
distribution narrower. However, to match cases in a simulated retrieval database with cloud side 
measurements we need to keep the number of parameters describing the relevant information about 
3D cloud structure as few as possible. As the next step, different combinations of radiances in our 
simulated retrieval database will be tested. 

Obviously, the retrieved values of droplet effective radius will correspond to droplets 
located not far (less than 1 km) from the cloud’s outer walls. However, in situ observations 
showed that the effective radius may remain constant for any given level in the cloud. In theses 
cases, retrieving effective radius near the cloud edges will give us information of the microphysics 
occurring in the cloud’s core. These features are discussed in more details by Mart ins et al. 
(2006). 

Acbrowledgments. This work was supported by NASA Goddard Space Flight Center New 
Opportunities Office. We thank Anthony Davis, Lorraine Remer, Steve Platnick, and Juan Roman 
for stimulating discussions 


14 



References 


Andreae M.O., D. Rosenfeld, P. Artaxo, A. A. Costa, G.P. Frank, K.M. Longo, and M.A.F Silva- 
Dias, 2004: Smoking rain clouds over the Amazon, Science, 303 (5662), 1337-1342. 

Blyth, A. M., and J. Latham, 2001: A climatological parameterization for cumulus clouds. J. 
Atmos. Sci., 48, 2367-2371. 

Breon F. M., and P. Goloub, 1998: Cloud droplet effective radius from spaceborne polarization 
measurements. Geophys. Res. Let., 25, 1879-1882. 

Breon F. M., and M. Doutriaux-Boucher, 2005: A comparison of cloud droplet radii measured 
from space. IEEE Trans. Geosci. Remote Sens., 43 (8), 1796-1805. 

Chang, F.-L. and Z. Li, 2002: Estimating the vertical variation of cloud droplet effective radius 
using multispectral near-infrared satellite measurements. J. Geophys. Res., 107, D15, 
doi: 1 0, 1 029/200 1 JD000766. 

Davies, R., 1978: The effect of finite geometry on the three-dimensional transfer of solar 

irradiance in clouds. J. Atmos. Sci., 35, 1712-1725. 

Davies, R., 1984: Reflected solar radiances from broken cloud scenes and the interpretation of 
scanner measurements. J. Geophys. Res., 89, 1259-1266. 

Davis, A. B., 2002: Cloud remote sensing with sideways-looks: Theory and first results using 
Multispectral Thermal Imager (MTI) data, in S.P.I.E. Proceedings, v. 4725, Eds. S.S. Shen 
and P.E. Lewis, S.P.I.E. Publications, Bellingham (Wa). 

Evans, K. F., 1998: The Spherical Harmonics Discrete Ordinate Method for Three-Dimensional 
Atmospheric Radiative Transfer, J. Atmos. Sci., 55, 429-446. 

Evans, K. F., S. J. Walter, A. J. Heymsfield and G. M. McFarquhar, 2002: The Submillimeter- 
wave cloud ice radiometer: Simulations of retrieval algorithm performance. J. Geophys. Res., 
107, D. 3, 10. 1029/200 1JD000709. 

Evans K.F. and A. Marshak, 2005. Numerical Methods in Three-Dimensional Radiative Transfer. 
In: A. Marshak and A.B. Davis, [Eds], “ Three-Dimensional Radiative Transfer in Cloudy 
Atmospheres” , Springer, 243-282. 

Iwabuchi H, and T. Hayasaka, 2002: Effects of cloud horizontal inhomogeneity on the optical 
thickness retrieved from moderate-resolution satellite data. J. Atmos. Sci., 59, 2227-2242. 

McFarlane, S. A., K. F. Evans, and A. S. Ackerman, 2002: A Bayesian algorithm for the retrieval 
of liquid water properties from microwave radiometer and millimiter radar data. J. Geophys. 
Res., 107, D16, 10. 1029/2001 JD001011. 

Marchuk, G., G. Mikhailov, M. Nazaraliev, R. Darbinjan, B. Kargin, and B. Elepov, 1980: The 
Monte Carlo Methods in Atmospheric Optics, 208 pp., Springer-Verlag, New-York (NY). 


15 



Marshak, A., A. Davis, R. F. Cahalan, and W. J. Wiscombe, 1994: Bounded cascade models as 
non-stationary multifractals. Phys. Rev. E, 49 , 55-69. 

Marshak, A., A. Davis, W. J. Wiscombe, W. Ridgway, and R. F. Cahalan, 1998: Biases in 
shortwave column absorption in the presence of fractal clouds. J. Climate , 11 , 431-446. 

Marshak, A. and A.B. Davis, [Eds], “ Three-Dimensional Radiative Transfer in Cloudy 
Atmospheres ”, Springer, 686 p. 

Martins J. V., Y. Kaufman, A. Marshak, D. Rosenfeld, L. Remer, and I. Koren 2006: Remote 
sensing of the aerosol effect on the vertical profile of cloud microphysics and 
thermodynamics, J. Geophys. Res., (to be submitted). 

McKee, T. B., and S. K. Cox, 1974: Scattering of visible radiation by finite clouds. J. Atmos. Sci., 
31, 1885-1892. 

Nakajima, T. Y., and M. D. King, 1990: Determination of the optical thickness and effective 
particle radius of clouds from reflected solar radiation measurements - Part I, Theory. J. 
Atmos. Sci., 47 , 1878-1893. 

Platnick S., 2000: Vertical photon transport in cloud remote sensing problems. J. Geophys. Res., 
105,22919-22935. 

Platnick, S., M. D. King, S. A. Ackerman, W. P. Menzel, B. A. Baum, J. C. Riedi, and R. A. Frey, 
2003: The MODIS cloud products: Algorithms and examples from Terra. IEEE Trans. 
Geosci. Remote Sens., 41 (2), 459-473. 

Prigarin, S., and A. Marshak, 2005. Numerical model of broken clouds adapted to observations. 
Atmospheric and Oceanic Optics, 18 , 256-263. 

Ramanathan, V., P. J. Crutzen, J. T. Kiehl, and D. Rosenfeld, 2001: Aerosols, climate and the 
hydrological cycle. Science, 294 , 2119-2124 

Rosenfeld, D., 2000: Suppression of rain and snow by urban and industrial air pollution, Science, 
287 , 1793-1796. 

Rosenfeld, D., and C. W. Ulbrich, 2000: Deep convective clouds with sustained supercooled liquid 
water down to -37.5°C, Nature, 405 , 440-442. 

Schertzer, D., and S. Lovejoy, 1987: Physical modeling and analysis of rain and clouds by 

anisotropic scaling multiplicative processes. J. Geophys. Res., 92 , 9693-9714. 

Stamnes, K., S.-C. Tsay, W. J. Wiscombe and K. Jayaweera, 1988: Numerically stable algorithm 
for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered 
media. Appl. Opt., 27 , 2502-2512. 

Stephens, G. L., 1994: Remote Sensing of the Lower Atmosphere. An Introduction, Oxford 
University Press, pp. 523. 


16 



Stephens, G. L., D. Vane, R. Boain, G. Mace, K. Sassen, Z. Wang, A. Illingworth, E. O’Connor, 
W. Rossow, S. Durden, S. Miller, R. Austin, A. Benedetti, C. Mitrescu and CloudSat Science 
Team (2002). The CloudSat mission and the A-Train: A new dimension of space-based 
observations of clouds and precipitation. Bull. Amer. Metereol. Soc., 83, 1771-1790. 

Varnai, T., and Marshak, A., 2002. Observations and analysis of three-dimensional radiative 
effects that influence MODIS cloud optical thickness retrievals. J. Atmos. Sci., 59, 1607- 
1618. 


17 



Figure captions 


Figure 1. Reflectance from a single cloud at two wavelengths: 0.67 pm (solid symbols) and 2.1 pm (empty 
symbols). Cloud height h = 2 km, cloud width L - 12.8 km, droplet effective radius, r t = 10 |im, SZA 9 0 =60° (a) A 
schematic illustration of illumination and viewing angles. Negative x correspond to reflectances from ‘cloud side’ while 
positive x correspond to reflectances from ‘cloud top’; (b) 0 = 60°; cloud optical thickness x= 160, 80, 40 and 20; (c) 
x= 80, 0=70°, 60°, and 45°. 

Figure 2. Reflectance from a single cloud with a variable droplet effective radius. Cloud height h = 4 km, cloud size 
L = 6.5 km, flat cloud top, t= 80, 0 O = 60°, 0= 45°. Droplet effective radius r e increases linearly with height from 5 
to 25 pm. (a) 0.67 urn; (b) 2.1 pm. Cloud edge is indicated by arrow atx = 7.5 km. Reflectance from cloud top is at 
the right side from the cloud edge while reflectance from cloud side is at the left. Dots indicate ‘measurements’ at ,v = 
0.25 km resolution. 

Figure 3. A realization of cloud stochastic model that has a given power-spectral exponent, mean, and standard 
deviation, (a) optical depth filed; (b) cloud top height field. 

Figure 4. An example of cloud microphysics, (a) liquid water content, LWC; (b) number of drops, N ; (c) droplet 
effective radius, r e ; (d) extinction coefficient, a esl . 

Figure 5. Reflectance from a realization of a stochastic cloud model with constant cloud optical (x=80) and geometrical 
thicknesses (h= 4 km). Left column: 0.67 pm; right column: 2.1 pm; 0 o =6O°; 0=0°, 20°, 45°, and 70°, from top to 
bottom. Note different color scales for left and right columns. 

Figure 6. A scatter plot of 2.1 pm reflectances vs. 0.67 pm reflectances based on 20 cloud fields generated by 
stochastic cloud model described in section 3. Parameters of the model are the following: mean cloud optical thickness 
= 80, mean cloud height = 4 km, spectral exponent = 2.0, standard deviations = 16 for the optical thickness and 1 km 
for the cloud height, cloud fraction = 0.5, 0 o =6O°, 0=70°. Particles smaller than 15 pm are water droplets while 
particles larger than 15 pm are ice. (a) Reflectances from both cloud sides and cloud tops. Dash lines indicate fixed 
0.67 pm reflectances (±0.03) used in Fig. 7. (b) Reflectance from cloud sides, (c) Reflectance from cloud tops. 

Figure 7. Histograms (number of cases vs. effective radius) obtained from Fig. 6. Values of reflectance at 0.67 pm 
were set to 0.83, 1.03, 1.22, and 1.50 with a ±0.03 window, on panels (a), (b) (c) and (d) respectively. Values of 
reflectance at 2.1 pm have a ±0.02 window. 

Figure 8. Histograms of reflectances at 0.67 pm and 2.1 pm conditional the effective radius equal to (a) 6-7 pm and 
(b) 14-15 pm. Plot is based on 20 realizations of the stochastic cloud model described in section 3. Parameters of the 
model are the same as in Fig. 6. 
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Figures 



Figure 1. Reflectance from a single cloud at two wavelengths: 0.67 pm (solid symbols) and 2.1 pm (empty 
symbols). Cloud height h = 2 km, cloud width L = 12.8 km, droplet effective radius, r t = 10 um, SZA 9 0 =60° (a) A 
schematic illustration of illumination and viewing angles. Negative x correspond to reflectances from ‘cloud side’ while 
positive x correspond to reflectances from ‘cloud top’; (b) 0 = 60°; cloud optical thickness %- 160, 80, 40 and 20; (c) 

T= 80, 6 = 70°, 60°, and 45°. 
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Figure 2. Reflectance from a single cloud with a variable droplet effective radius. Cloud height h = 4 km, cloud size 
L = 6.5 km, flat cloud top, %= 80, 0 O = 60°, 0= 45°. Droplet effective radius r e increases linearly with height from 5 
to 25 pm. (a) 0.67 pm; (b) 2.1 pm. Cloud edge is indicated by arrow at x = 7.5 km. Reflectance from cloud top is at 
the right side from the cloud edge while reflectance from cloud side is at the left. Dots indicate ‘measurements’ at 5 = 
0.25 km resolution. 
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Figure 4. An example of cloud microphysics, (a) liquid water content, LWC; (b) number of drops, N\ 
effective radius, r e ; (d) extinction coefficient, a cxt . 
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Figure 5. Reflectance from a realization of a stochastic cloud model with constant cloud optical (t= 80) and geometrical 
thicknesses (/i=4 km). Left column: 0.67 pm; right column: 2.1 pm; 0 o =6O°; 0=0°, 20°, 45°, and 70°, from top to 
bottom. Note different color scales for left and right columns. 
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Cloud tops 


Figure 6. A scatter plot of 2.1 pm reflectances vs. 0.67 pm reflectances based on 20 cloud fields generated by 
stochastic cloud model described in section 3. Parameters of the model are the following: mean cloud optical thickness 
= 80, mean cloud height = 4 km, spectral exponent = 2.0, standard deviations = 16 for the optical thickness and 1 km 
for the cloud height, cloud fraction = 0.5, 0 o =6O°, 0=70°. Particles smaller than 15 pm are water droplets while 
particles larger than 15 pm are ice. (a) Reflectances from both cloud sides and cloud tops. Dash lines indicate fixed 
0.67 pm reflectances (±0.03) used in Fig. 7. (b) Reflectance from cloud sides, (c) Reflectance from cloud tops. 
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Figure 7. Histograms (number of cases vs. effective radius) obtained from Fig. 6. Values of reflectance at 0.67 (im 
were set to 0.83, 1.03, 1.22, and 1.50 with a ±0.03 window, on panels (a), (b) (c) and (d) respectively. Values of 
reflectance at 2.1 urn have a ±0.02 window. 







Figure 8. Histograms of reflectances at 0.67 ftm and 2.1 jam conditional the effective radius equal to (a) 6-7 iim and 
(b) 14-15 !im. Plot is based on 20 realizations of the stochastic cloud model described in section 3. Parameters of the 
model are the same as in Fig. 6. 





